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Abstract. Ground state of the two-dimensional hard-core-boson system in the presence of the quenched 
random chemical potential is investigated by means of the exact-diagonalization method for the system 
sizes up to L — 5. The criticality and the DC conductivity at the superfluid-localization transition have 
been controversial so far. We estimate, with the finite-size scaling analysis, the correlation-length and the 
dynamical critical exponents as v = 2.3 ± 0.6 and z = 2, respectively. The AC conductivity is computed 
with the Gagliano-Balseiro formula, with which the resolvent (dynamical response function) is expressed 
in terms of the continued-fraction form consisted of Lanczos tri-diagonal elements. Thereby, we estimate 
the universal DC conductivity as a c (uj — > 0) = 0.135 ± 0.01((2e) 2 //i). 

PACS. 75.10.Jm quantized spin models - 75.10.Nr spin glass and other random models - 75.40.Mg 
numerical simulation studies 

1 Introduction instance, the integer quantum hall effect is described in 

terms of successive metal-insulator (delocalization-localization) 
The scaling argument of Abrahams, Anderson, Licciardello transitions of dirty two-dimensional electron system with 
and Ramakrishnan § states that in two dimensions, in- the external magnetic field varied. In the above-mentioned 
finitesimal amount of quenched randomness should drive scalin S theor y' the random perturbation is appeared to be 
itinerant extended states to localize. That is, at the abso- marginal so that some unexpected factors, namely, the 
lute zero temperature, the conductivity should be vanish- magnetic field and the many-body interaction, would pos- 
ing, if there exist any randomnesses. These are, however, slD ly change the scenario, 
some exceptions where the above description fails. For 
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For instance, suppose that there exists an attractive nomcnon occurs at the ground state. This criticality is 
interaction among the electrons. At the ground state, the different essentially from that of the finite-temperature 
electrons would be unstable against the bose condensation transitions. In the path-integral picture for the partition 
so that the system would be in the superconducting phase, function, d-dimensional quantum system is regarded as a 
The disorder-driven localization from the superconducting (d+l)-dimensional classical system. The system-size along 
phase is apparently out of the scope of the conventional the imaginal-time direction is given by the inverse tem- 
localization theory, and has been studied extensively so perature, which is diverging at the ground state. Hence, 
far. In experiments, the transition is observed for metal- the critical fluctuation along the imaginary-time direction 
lie films [0 ,|L[|jf| , high-T c films Q and Josephson-junction contributes to the universality class significantly. This ex- 
arrays @,|. One of the main concerns is the conductivity tra contribution is characterized by the dynamical critical 
at the critical point: The experiments show that irrespec- exponent z, which is explained in the next section in de- 
tive of the samples examined, the localization transition tail. 

occurs at a universal condition. Namely, at the transition, m. ,• , ifiTjiriii i ,1 , .-i j 

' I he scaling argument ||12|,I14I| shows that the conduc- 

the electrical conductivity is found to be ~ (2e) 2 /h. (The . . .. -a ^ * 1 v , 

J \ / / \ tivity remains finite at tne onset of tne localization tran- 

parameter e denotes the charge of a single electron, and h ... . ,-, , , , 

r sition even at tne ground state. Moreover, tne argument 

denotes the Planck constant.) c , . . , , , 

' tor general dimensions yields tne dynamical critical expo- 
Localization transition of the absorbed helium by a nent being equal to the spacial dimension; z = d. We ex- 
porous media is considered p!l|,^2| to belong pre- plain the argument in the subsection 2.2. This prediction 



cisely to the same universality class mentioned above. (That is astonishing in the sense that there is no upper critical 

is, the critical exponents are identical. The electrical con- dimension in this critical phenomenon. (The prediction 

ductivity, however, does not make any sense for the latter, is confirmed rigorously for d = 1, where the bosoniza- 

because the helium atom is not charged.) The equivalence tion technique is available JT3| , p^[ . For numerical study of 

is based on the belief that the Cooper pair is formed al- the d = 1 criticality, see the article jl6| and references 

ready in the localization phase as well as in the supercon- therein.) Because of the absence of the upper critical di- 

ducting phase, and the essence of the transition is con- mension, the e-expansion scheme cannot be formulated, 

cerned only in the boson degrees of freedom and the site- Hence, for d = 2, numerical-simulation studies have been 

random chemical potential. This picture was found to be playing a crucial role so far. 

valid in one dimension |l3| . Runge |i~7j employed the exact-diagonalization method 

From a theoretical viewpoint, the criticality itself is to treat system sizes up to L = 4; the Hamiltonian is 

a matter of interest. It is notable that the present phe- given by eq. (0) shown afterwards. (When treating L = 5, 
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he reduced the particle density to n = 0.2.) He obtained random chemical potential. The Hamiltonian is given by, 

the estimates v = 1.4 ± 0.3, z — 1.95 ± 0.25 and a c = „. J v-^ ( + \ \ , T ,„ + 

H = --J^(a\a j + h.c.) +2fr i (2ojo i -l), (1) 

0.17 ± 0.01((2e) 2 //i), and claimed that finite -size correc- w ' 

. „ .. , . , T „ . ,, ,, where the operators ia^ajl obey the hard-core boson 
tion prevents definite conclusion. Monte-Carlo metfiod fias 1 ,J 

been used in many studies. Makivic, Trivedi and Ullah jl8| s ^ a ^ ls lcs ' 

obtained the estimates v = 2.2 ± 0.2, z = 0.5 ± 0.05 and [a it aj] — [a\,a^] — [a i} a]] =0 (i ^ j), 

a c = 1.2±0.2. Batrouni, Larson, Scalettar, Tobochnik and { ai ,at} = l and a,a t = a\a\ = 0, (2) 

Wang [|l9) obtained the estimate a c = 0.45 ±0.07. Wallin, 

and denotes the summation over all nearest neigh- 

S0rensen, Girvin and Young [^0| obtained v = 0.9 ±0.1, 

bors. The site-random chemical potentials {Hi} distribute 

z = 2.0±0.1 and a c = 0.14±0.03. S0rensen, Wallin, Girvin 

uniformly over the range [— \/3A\; namely, the mean 
and Young (2| obtained v = 1.0 ± 0.1, z = 2.0 ± 0.1 and s 

deviation is given by ^[Hf]^ = A. The particle density 

a = 0.14 ± 0.03. Zhang, Kawashima, Carlson and Guber- 

n = N/L 2 is fixed to be one half throughout this paper. 

natis jH obtained v = 0.9 ± 0.1 and z = 2.0 ± 0.4. The 

As is mentioned above, the model (Q) is believed to 

estimates are still rather scattering. In particular, comput- 

describe the physics of the superconductivity-localization 

ing the conductivity more reliably would be of practical 

transition as well as the superfluid-localization transition. 

importance, because it can be compared with the exper- 

This belief is based on the picture jOlpJI that the boson 

imentally observed value <r c ~ (2e) 2 /h mentioned above. 

describes the Cooper pair, so that the boson charge should 

The Monte-Carlo method enables one to treat larger sys- 

be the twice of the single electron charge; e* = 2e. The 

terns than those with the exact-diagonalization method. 

model ([j]) is investigated very extensively in the above 

In the quantum Monte Carlo simulation, however, because 

mentioned article Jl^j . Here, we attempt to improve this 

the dynamical critical exponent is predicted to be z — 2, 

work through utilizing some recent developments, the al- 

in order to simulate ground-state property, the imaginary- 

gorithm for computing the response function (resolvent) 

time system size should be enlarged rapidly; namely, it 

p3| , the estimate scheme of the dynamical critical expo- 
should be kept quadratic in terms of the real-space sys- 

ncnt [ p4[ , and parallel supercomputers which enable one 

tern size at least. 

to treat larger system (L = 5). These improvements are 
described in respective subsections of Section |^. 

The rest of the paper is organized as follows: In the 
In the present paper, we simulate the two-dimensional next section, we review some notions relevant to the present 
bose system on the L x L square lattice under the peri- study: First, we show the viewpoint from which the model 
odic boundary condition in the presence of the quenched (0) is regarded as the quantum XY model |25|| . Nature of 
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the superfluid-localization transition is interpreted in the in terms of the spin operators, 

language of the quantum spin system. Then, we summa- , — , , — , 

H= -jJ2( S i S j+ S ! S j)+ 2 J2 HiS i- (3) 

rize the scaling argument |l2|,[l4| describing the critical- (v) i 

mi • u t c i i-i Now, perspectives developed for the quantum spin sys- 
lty. Ine argument yields some scaling formulae wmcn are 

r , . ,. -ii,t;i tern become available. Kishi and Kubo [E6l showed rig- 
useful m analyzing our finite-size numerical data. In tne Lt - J 

pn , orously that (without the random magnetic field) in the 

section pi our numerical simulation results are presented. j \ o / 

,, T , j ;i i n o i n r- j thermodynamic limit, the ground-state magnetism is spon- 

We estimated the critical exponents as v = 2.3 ± 0.6 and J 



z = 2. At the critical point, we estimate the conductivity 
as er c = O.f 35 ± O.Of ((2e) 2 /h). These new results are sum- 



taneously broken; 

(m 2 Yy ) * 0, (4) 



marized in the last section in comparison with previous where mxy = XwC^f > ^ ^ s expected that for suf- 
results. ficiently strong random fields, the long-range magnetism 



2 Review — equivalence to the XY model 
and scaling argument 

In this section, we summarize several important aspects 
about the model (Q). First, we introduce the equivalence 
between the model ([j]) and the quantum XY spin system 
p5f . This equivalence reveals nature of both the superfluid 
and the localization phases in the language of the spin sys- 



would be disturbed, (m 2 XY ) = 0- Hence, we see that a 
phase transition between the XY and the random-field 
phases would exist at a certain random-field strength. 

What does the existence of the XY order (Q) stand for 
in the boson language? We show that it stands for the bose 
condensation (superfluidity) : If the number of the particles 
condensing at the k = state is of the system-size order, 
the state is regarded as being of superfluid. That is, the 
quantity, 



k=0 



tern, and provides intuitive picture of the phase transition. |^supcr| 2 = -j-p, ^ ( L & ^ ^ j ( L 6 * ^ j 

Finally, we review the scaling argument 14 , whose for- (5) 

mulae are used in the analyses of our numerical data in works as an order parameter of the superfluidity. The or- 

Section ^. der parameter is expressed in term of the spin language, 

Itf'supcrl 2 = ( m xi-)- As is explained above, this remaines 

„..„.. . , »„, ■ . . i , finite. And so, the superfluidity develops actually at the 

2.1 Mapping to the quantum XY model with the J 

, . ,. .. ground state of the Hamiltonian (TO) at A = 0. It is quite 

random magnetic field ° \o» ^ 

natural that the gauge degree of freedom, which is spon- 
Because the S = 1/2 ladder operators {Sf , Sf} obey the taneously broken in the superfluid phase, is related to the 
same algebra eqs. (||) p5| , the Hamiltonian (|l|) is expressed in-plain spin rotator. As is introduced in Section [|, the 



Yoshihiro Nishiyama: Two-dimensional ground-state superfluid-localization transition 5 

essence of the superfluid-localization transition is believed general, for quantum random systems, the real-space cor- 

to be concerned only in the in-plain rotator degrees of relation develops less robustly than the imaginary-time 

freedom, and the site-random perturbations conjugate to one does. This anisotropy is one of the significant char- 

them. acteristics of the random quantum critical phenomena, 

resulting in various new exotic universality classes. The 
2.2 Scaling argument anisotropy is characterized by the dynamical critical ex- 
Here, we introduce a scaling argument [ ^2|Jl^1 which dc- P onen t z - 

,, r,, a -j i v • The superfluid density p s — the spin stiffness in the 
scribes the criticahty of the superfluid-localization transi- 

, . „. . , , , „ , „ , language of the spin system — is defined as the clastic 
tion. Ine argument yields various useful formulae to an- 

, i j i rni ir i j constant in terms of the real-space gauge twist [E7|, 
alyze numerical data, lhe argument itself, however, does 

not yield any conclusions to estimate the critical expo- „ Ps ,^ „o / 7 \ 
ncnts quantitatively. Any analytical analyses to estimate 

, . , . . „ , • m i • • i ■ Note that the quantity works as an order parameter of 
the critical exponents face a difficulty as is mentioned m 

TJ t^;i r;- ;• ,, the superfluidity. It is furthermore expressed in the form, 
Introduction, .tor tne purpose of estimating tne exponents 

ii. i i i / ~ % (y-) 2 , where denotes the total gauge twist. Using 

quantitatively, numerical simulation lias been playing a, z y ^ ' 

. , , the form (o), we obtain, 

crucial role. ^ 

The scaling hypothesis states that (the singular part ( T 2 

Pa ~ J " 

of) the free energy per unit volume for the system size (d+z 2) ~ ( ^ P\ 

= ^ ^ ' ir ) ^ 

L and the inverse temperature j3 should be given in the _ ^ ^ 

form, \Cr £,r ) 

f(L,{3) ~ f ! (^) This scaling form is used in our numerical-data analysis. 

where £ r and £ T denote the real-space correlation length The compressibility n is defined, on the other hand, as 

and the imaginary-time one, respectively. (Note that in the elastic constant of the imaginary-time gauge twist, 

the path-integral viewpoint, the partition function of a d- f ~ f (<9 r Q) 2 . Through the similar arguments as the above, 

dimensional quantum system is regarded as that of a (d + we obtain the scaling formula for the compressibility, 

l)-dimensional counterpart.) These are expressed in term ^ f L (3 



K~^ d - Z) f[f,f )■ (9) 
of the deviation from the critical point 5; £ r ~ 8 v and 

£ T ~ £*. The parameter v denotes the correlation- length Fisher, Weichman, Grinstein and Fisher predicted that 
critical exponent, and z denotes the dynamical critical ex- the formula, 

ponent. These are estimated in Section ^, numerically. In z = d, (10) 
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should hold for any dimensions: The compressibility is fi- bers N = [L 2 /2] and [L 2 /2] — 1; the bracket [• • •] denotes 

nite in both phases beside the transition point. Hence, it the Gauss notation. The data for n = 0.5 are obtained 

would be kept to be of the order unity even at the crit- through interpolating these two sets of data with use of 

ical point as well, so that we obtain the equality (|l0|). the relation Q(n) = a(n — 0.5) 2 + c; physics is symmetric 

In Section |^, we confirm their prediction with use of the in terms of n = 0.5 (particle-hole symmetry). 
Rieger- Young method |24j . 

Finally, we explain how the above scaling argument 3.1 Criticality of the superfluid-localization transition 



concludes that the conductivity would be finite at the crit- 
ical point. The scaling argument juj yields the following 
formula for the AC conductivity, 



Here, we determine the location of the superfluid-localization 
transition point and the correlation-length critical expo- 
nent. We use the language of the spin system, which is 
o'M ~ £ (d+Z 2) Ps(^£r)/w- (11) explained in the subsection O. That is, from this view- 



point, the transition is characterized by the disappearance 
of the XY (in-plain) magnetic order. 



Assuming d — 2, and the scaling function behaves as 

p s (x) oc x in the vicinity of the critical point x — > oo, we 

, , . ,1 • i •;• In Fig. hi we plotted the square of the in- plane magne- 

sec that the conductivity does not have any singularities at a 

,, . . , . i j.i r n -± rn tization m? = [(mvv)lav against the randomness, where 

tne critical point, and therefore remains finite, lo summa- aih o 

,i n . , ., ., the bracket [• • -lav denotes the random-sampling average, 

rize, the scaling argument states that there is a possibil- L J ' 

i , . n , a ■ , and (• • •) denotes the ground-state expectation value. The 

lty that tne conductivity remains finite at tne supcrnuid- 

, ,. , . ... , . , ... , . , , . „ random-sample numbers are 2048, 2048 2048 and 608 for 

localization critical point, lne quantitative evaluation ot 



er c , however, lies out of the scope of the argument. Here 

we compute the conductivity with use of the Gagliano- 

n . i .cm • ,, , . Fig. 0, we plotted the Binder parameter pal of the in-plane 

Balseiro method 23 m the next section. 13 



L = 2, 3, 4 and 5, respectively. We see that the in-plane 
magnetization becomes suppressed by the randomness. In 



magnetic order, 

3 Numerical results U = l- ^HTdT ■ ( 12 ) 



In this section, we present our numerical results. In or- The Binder parameter is invariant with respect to the 

der to diagonalize the Hamiltonian ([!]), we employed the system sizes at the critical point. It is enhanced (sup- 

Lanczos method. We have fixed the particle density to pressed) in the order (disorder) region as the systems size 

be one half n(= N/L 2 ) = 0.5, and treated the system is enlarged. In Fig. |2| We observe an intersection point at 

sizes up to L — 5. For those system sizes of odd L, we A w 1.2. Namely, in the region A < 1.2, the superfluidity 

proceeded the sets of simulations for the particle num- persists against the random chemical potential, whereas 
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Fig. 1. Square of the in- plain magnetization m 2 is plotted 
against the randomness A. The random-sample numbers are 
2048, 2048, 2048 and 608 for L = 2, 3, 4 and 5, respectively. 
The XY magnetic order becomes suppressed as the random 
magnetic field is strengthened. 



0.5 - 



0.4 



0.3 
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Fig. 2. Binder parameter ( p^ ) of the in-plain magnetization is 
plotted against the randomness. The intersection point of the 
curves indicates the location of the transition point. 



in A > 1.2, the particles arc localized, and the long-range 
gauge coherence is lost. 

In order to estimate the transition point and the correlation- 
length exponent precisely, we analyze the above data by 
means of the finite-size-scaling theory According to the 
theory, the Binder parameter (dimensionless quantity) obeys 
the form U = U ((A - Z\ c )i 1 / l/ ) in the vicinity of the crit- 
ical point A c . Namely, the data (A — A^L 1 ^-!/ should 
collapse along a universal curve irrespective of the system 
sizes. In other words, the critical point A c and the expo- 
nent v are adjusted so that the scaled data could form a 
universal curve. The degree to what extent the data col- 
lapse is measured by the Kawashima-Ito "local-linearity 
function" J29j which is explained in Appendix. In Fig. |^, 
we show the scaling plot for the data shown in Fig. 0. In 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



L=2 
L=3 
L=4 
L=5 



-1 -0.5 0.5 

(A-1.25)L 1/2 '° 



Fig. 3. Scaling plot for the data shown in Fig. g. The scaling 
analysis yields the best-fit estimates A c = 1.25 and v — 2.0. 



consequence, we obtain the estimates A c — 1.25 ±0.1 and 
v = 2.3±0.6. In order to see whether the correction to the 
finite-size scaling exists, we analyzed the data for L = 3, 
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4 and 5 similarly; omitted the data for L — 2. This anal- 
ysis yields the best-fit estimates A c = 1.27 and v — 2.3. 
Therefore, we conclude that there is little correction to the 
finite-size scaling; in other words, the system sizes treated 
here reach the scaling region. Our estimate v = 2.3 ±0.6 is 
somewhat different from that of ref. |l7| (y — 1.4±0.3) de- 
spite of the fact that we considered the same model also 
using the exact-diagonalization technique. This discrep- 
ancy may originate in the criteria utilized to appreciate 
the scaling-data collapse, and is discussed in the next sec- 
tion. 

According to the formula (||), if we assume that the 
dynamical critical exponent is equal to the spacial dimen- 
sion, namely, z = d = 2, the scaled spin stiffness L 2 p s 
should be invariant at the critical point. The spin stiffness 
is given by, 

d 2 E g (0) 



dO 2 



0=0 



(13) 



where the angle denotes the boundary gauge twist a\a,j — > 
e l0 a\aj, and E g denotes the ground-state energy. In Fig. 
f|, we plotted the scaled spin stiffness L 2 p s . In fact, we 
observe at A « 1.3, all the curves intersect. In conse- 
quence, founded on the scaling formula (||), the so-called 
"generalized Josephson relation," we see that the dynam- 
ical exponent would be equal to two. The same reasoning 
founded on this relation has been reported in the liter- 
atures 17,1^,^;^2|. In this paper, in Section 3.3, we es- 
timate the dynamical critical exponent with use of the 
Rieger- Young method, which would be more straight- 
forward. 



-J 




Fig. 4. Scaled spin stiffness L 2 p B is plotted. Because the curves 
intersect at a point, the prediction z = d is appeared to hold. 

3.2 Electrical conductivity 

We evaluate the dynamical conductivity at the critical 
point A c = 1.25 estimated in the above subsection. The 
DC conductivity is conjectured to be universal at the crit- 
ical point as is explained in Introduction. The dynamical 
conductivity is given by the current-current time correla- 
tion, 



<j{ijj) = Re 



1 1 f°° 

J^T 2 1 dteIwt <[ J -W> J *]> 



(14) 



where the current is given by J x — s S x Oj +s a,j. 

In the previous work ]l7j |, the conductivity is calculated 
through Fourier-transforming the current-current correla- 
tion numerically, whose time correlation is evaluated with 
use of the Suzuki- Trotter-decomposition approximation. 
We evaluate the resolvent form of eq. (E3) directly, 



<t(lj) = Re 



/ i 



\hwL 2 

h 
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Now, we are free from the discrete numerical integration 
error. (The conductivity is a linear combination of the 
delta- function peaks as is apparent from eq. ([l5]), which 
might be suffered significantly from the numerical dis- 
cretization error.) Some might wonder that the inverse 
matrix of the total Hamiltonian in eq. ([l5]) cannot be com- 
puted; this is true. The expectation value of the inverse of 
the Hamiltonian is, however, evaluated with use of the 
Gagliano-Balseiro continued- fraction formula p3[ , 



/o 



1 



z-H 



fo 



(/o|/o) 



z - a Q 



1L 



where the coefficients are given by the Lanczos tri-diagonal 
elements, 



l/i+i) =H\f t )-a l \f l )-f3f\f^ 1 ), (17) 

<n = (fi\n\fi)/{fi\fi), 

01 = (fi\fi)/(fi-i\fi-i) Wa = 0). 

We have evaluated the dynamical conductivity (|l5|) by 
means of this formula. 

In Fig. we show the conductivity for a certain ran- 
dom sample with A = 1.25 and L = 5. We see that the 
conductivity consists of delta-function peaks. The delta 
function peak is broadened into the Lorentz form through 
the substitution to — > u — 0.2L The concept of the dissi- 
pation becomes subtle for finite-size system. For instance, 
the DC conductivity is vanishing. This is due to the pres- 
ence of the finite-size energy gap above the ground state. 
Only in the thermodynamic limit, The conductivity in the 
vicinity of the zero frequency w = would emerge. (In the 
quantum Monte Carlo simulation, the (real-)time corre- 




0.000 



(16) Fig. 5. Electrical conductivity is evaluated for a random sam- 
ple with A = 1.25 and L = 5. We see that the conductivity 
consists of delta-function peaks; the delta-function is broad- 
ened into the Lorentz form with the width r\ = 0.2. 



lation cannot be computed; the temperature (imaginary 
time) correlation function is computed instead. The latter 
is free from the delta- function singularities, because the 
poles exist only along the real-frequency axis. An analyt- 
ical continuation to the real frequency, however, should 
be performed in order to estimate the dynamical response 
function such as eq. (|l4|). This is extremely difficult.) 

In Fig. ||, we show the random-averaged conductivity 
for A = 1.25; the random-sample numbers are 16384, 2048 
and 960 for L ~ 3, 4 and 5, respectively. We see that 
the conductivity increases as the frequency is reduced. 
In the vicinity of the static point uj ~ 0, however, the 
conductivity drops rapidly due to the reason mentioned 
above (finite-size effect). (For comparison, we show the 
Fourier-transformed results jl7|; a « 0.055 for L = 3 and 
a w 0.105 for L — 4. These are a bit larger than ours.) 
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0.125 



0.1 00 




0.000 



Fig. 6. Electrical AC conductivity is evaluated with the 
Gagliano-Balseiro formula ( |l6| ) for A = 1.25; the random- 
sample numbers are 16384, 2048 and 960 for L — 3, 4 and 
5, respectively. 

We estimate the DC conductivity as the maximal value 
of the AC conductivity. The DC conductivity shows large 
system-size dependence. In Fig. 0, we depict the 1/L 2 ex- 
trapolation of the conductivity. (The power 2 is chosen for 



0.14 



0.12 



0.10 



o 0.08 - 

Q 

b 

0.06 - 



0.04 



0.025 0.05 0.075 

1/L 2 



0. 1 25 



Fig. 7. 1/L 2 extrapolation of the DC conductivity. 



the same reasoning as that in the paper |17|| .) The plots 



align. We stress that the present new data for L — 5 is 
crucial to confirm the validity of the 1/L 2 extrapolation. 
We obtained the extrapolated conductivity with the least- 
square fit as CT C = 0.135 ± 0.01 ((2e) 2 /^)- 

3.3 Dynamical critical exponent 

According to the conjecture introduced in Section \}.2\ 
the dynamical critical exponent is given by z = 2 at 
the superfluid-localization transition. This conjecture has 
been confirmed numerically jl7],^0|,[22j with the help of the 
generalized Josephson relation (0) for the superfluid den- 



sity. In fact, we demonstrated in Section 3.1 that the su- 
perfluid density is well described in terms of the gener- 
alized Josephson relation and the assumption z = 2. In 
this subsection, we utilize the Rieger- Young formula [ p4| 
for the first time in order to estimate z. We think that 
the scheme is suitable for the exact-diagonalization simu- 
lation, and gives z straightforwardly. 

In Fig. ||, we show the probability distribution of the 
first energy gap AE; A — 1.25, L — 5 and 960 random 
samples. According to the Rieger- Young argument, from 
the low-energy tail of the distribution, the dynamical crit- 
ical exponent is extracted: The probability of a certain 
energy gap may be proportional to the spatial volume 
P (x L d ■ AE X , where the exponent A describes the low- 
energy tail. On the other hand, the finite-size scaling the- 
ory states that any quantities should be expressed in the 
form P = P(L/£ T ) = P{L/£ /z ); see the subsection |]§ 
With use of the formula AE ~ 1 /£ T , the distribution turns 
out to be a function of L(AE) 1 / 2 . Through collating this 
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Fig. 8. Probability distribution of the first energy gap for the 
system with A — 1.25 and L = 5. 

fact with the above form, we obtain the relation A = d/z, 
which is the Rieger- Young relation. 

In Fig. ||, we see that the low-energy tail is, in fact, 
well described by the exponent A = 1; namely, we obtain 
the estimate z = d = 2. 

4 Summary and discussions 

We have investigated the two-dimensional hard-core bo- 
son system in the presence of the random chemical poten- 
tial, whose Hamiltonian is given by (Q). We diagonalized 
the Hamiltonian numerically for the system sizes up to 
L = 5. The present new data for L = 5 show no unex- 
pected (irregular) behavior in the finite-size-scaling anal- 
ysis. That is, the system sizes treated here reach scaling 
regime. This was not so certain, and in fact worried about 
in the previous study. We observed that the disorder- 
driven superfluid-localization transition takes place at the 
critical randomness A c = 1.25 ± 0.1. We estimated the 



correlation-length exponent and the dynamical critical ex- 
ponent as v = 2.3 ± 0.6 and z = 2, respectively, with 
the finite-size scaling analysis. By means of the Gagliano- 
Balseiro method, we computed the dynamical conductiv- 
ity at the critical point. Thereby, we obtained the DC 
conductivity a c = 0.135 ± 0.01((2e) 2 /h). 

Our estimate of the critical point is accordant with 
that of Runge pj] : By means of the same numerical method 
as ours, Runge obtained A c ~ 1.15 through the analysis 
of p s , A c ~ 1.25 through the Binder parameter (the same 
as ours) and A c ~ 1.1-1.3 through (m^- y ). The present 
new simulation for larger system size (L = 5) confirmed 
that the system size is surely in the scaling regime. On 
the other hand, by means of the quantum Monte Carlo 
method it has been reported that the transition point 
locates at A c = 1.43 ± 0.06 @ and 0.72 Q (in our 
definition). These estimates are based on the generalized 
Josephson relation (^]) and an assumption of z; the former 
assumed z — 2, whereas the latter assumed z = 0.5. We 
think that the Binder parameter of the gauge order ([l2|), 
which is readily computed with the exact-diagonalization 
scheme, shows less correction to the finite-size scaling. 

Secondly, we discuss the correlation-length exponent 
v. It should be made clear why our estimate v = 2.3 ± 0.6 
differs from that of rcf. JTtJ [y = 1.4 ± 0.3) despite of 
the fact that we considered the same model by means of 
the same numerical technique except for our extension to 
L = 5. Crucial difference might be the criteria utilized 
to appreciate the data collapse in the scaling analysis. As 
is explained in Appendix, we used a quantitative scheme 
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which also gives the estimate of the error margin. The vital in order to confirm the validity of the 1/L 2 extrap- 

error margin is, however, determined through taking ac- olation. Quantum Monte Carlo method is less efficient in 

count of only the statistical error, but not the correction to computing dynamical quantity. (In the quantum Monte 

finite-size scaling. The latter is hardly appreciable quanti- Carlo simulation, the temperature (imaginary time) cor- 

tatively. Hence, there is a possibility that the error margin relation rather than the time correlation is evaluated in- 



would be somewhat larger. The estimate of ref. |17 lies a stead. The temperature correlation is fitted by an analyt- 

bit outside of our error bar, and, in fact, the scaling plot ical function so as to yield the time correlation through 

based on the assumption v = 1.4 does not show distinct the Wick rotation.) The present calculation with the help 

failure as far as one can see. We believe, however, that our of the Gagliano-Balseiro method compensates the disad- 

quantitative criterion would be objective and possibly less vantage. 

biased, although the present error margin might be rather In the paper [|2] , the authors claim that the limit 

large. At least, our data shown in Fig. || seem to exclude u) — > at T = 0, which is used in the present paper, fails 

the conclusion v ~ 1. In fact, Chayes et al. claimed an to give the DC conductivity relevant to the experimental 

inequality v>2/d, p(i| , |3l"[ which casts doubt on v < 1. transport measurements. According to their theory, the 

Next, we turn to discuss the dynamical critical ex- AC conductivity shows a Drude-like peak with the width 
ponent. We have used the Rieger-Young relation for the Hence, at T -> 0, the peak vanishes to a single 

first time to estimate the dynamical critical exponent. The P oint so as to make the AC conductivity singular (dis- 

analysis showed a clear evidence that z{= d) = 2 holds, continuous) at uj = 0. Through taking account of this, 

Because the previous estimates founded on the general- the authors speeded in explaining the experimentally 

ized Josephson relation (§) indicate the same conclusion observed critical conductivity « (2ef/h. Our numerical 

as well, we believe that the prediction z = d is established simulation is rather incapable of observing the features 

fairly definitely in two dimensions. the y P ro P osed - Thi « remained to be solved in future. 

Finally, we mention about the electrical conductivity. 

The present conclusion a c / ((2e) 2 /h) = 0.135±0.01 is com- Acknowledgement 
parable with the estimates 0.17 ±0.01 £L7| and 0.14 ±0.03 

iHTji TTr ,, ,,. . , ,. ,. ,. , , Our program is based on the subroutine package TIT- 

I P Ofl . We stress that the previous exact-diagonahzation data 

0. , „ . . . n . PACK ver. 2 coded by professor H. Nishimori. Numerical 

tor L = 3 and 4 are rather far from convergence, and 

,, r T . i i , i , , i i, r ,i simulations were performed on VPP/700 56 of the com- 

our result tor L — 5 reveled that the result of the system ' ' 

n , . puter center, Kyushu university, 

size is about to converge to a certain tnermodynamic-limit 

value. In fact, as is shown in Fig. 0, our new data L = 5 is 
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A Details of the present scaling analyses we use causes systematic error. Note that in the present 



We explain the details of our finite-size-scaling analyses, 
which we managed in Section [}] in order to estimate the 
transition point A c and the exponent v. We adjusted these 
scaling parameters so that the scaled data shown in Fig. 
^ form a universal curve irrespective of the system sizes. 
In order to see quantitatively to what extent these data 
align, we employ the "local linearity function" S defined 
by Kawashima and Ito p9| |: Suppose a set of the data 
points {(xi,yi)} with the error-bar {di(= 5yi)}, which we 
number so that Xi < Xi+\ may hold for i = 1, 2, • • • , n — 1. 
For this data set, the local-linearity function is defined as 



analysis, we do not have to assume any fitting functions. 

The scaling parameters are determined so as to mini- 
mize the local-linearity function. The error margin of the 
scaling parameter is hard to estimate. It is concerned with 
both the statistical error and the correction to the finite 
size scaling. The former error can be estimated through 
considering the statistical error of the function S. This 



function has a relative error of the order 1 / ^Jn — 2. (The 
local-linearity function is of the order (n— 2), and the error 



is given by \Jn ~ 2.) The correction to the finite-size scal- 
ing, on the contrary, is quite difficult to determine. Here, 
we consider only the statistical error in order to estimate 
S = ^2 w ( x i>y^ d i \x i - 1 ,y. l - 1 ,d l _ ll x i+ll y i+1 ,d t+1 ). the error marging of tne scamlg parameters. As the num- 



n-l 



(18) 

The quantity w (xj , yj , d 3 ■ | x l , y l , di , x k , yk , d k ) is given by 



, 2 

Vj - V 
A 



where 



and 



(x k - x 3 )yi - (xi - Xj)y k 



Xk Xi 



ber of the data points n is increased, the statistical error 
of S is reduced. The corrections to the finite-size scaling 
j^g) might increase instead. In the present analyses, we used 
twenty data in the vicinity of the transition point A c . An 
example of the plot S is shown in Fig. 0. We observe the 
minimum at v — 2.0. Taking into account of the statistical 
error, we estimated the critical exponent as v — 2.3 ± 0.6. 



(20) 
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